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^r*j ' Abstract 

o ■ 

■ The generating functional of the cyclic representation of the CCR (Canonical Com- 

mutation Relations) representation for the thermodynamic limit of the grand canon- 
ical ensemble of the free Bose gas with attractive boundary conditions is rigorously 
computed. We use it to study the condensate localization as a function of the homoth- 
ety point for the thermodynamic limit using a sequence of growing convex containers. 
The Kac function is explicitly obtained proving non-equivalence of ensembles in the 
condensate region in spite of the condensate density being zero locally 
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1 Introduction 



The interest in the phenomenon of standard Bose Einstein Condensation (BEC) revived in 
recent years due to the spectacular experimental work on Bosons in traps. We refer, e.g., 
to [9] and [12] for experimental and theoretical state of affairs. A renewed interest in old 
problems connected with the phase transition accompanying BEC is at order. The generic 
model for BEC is the free Bose gas as already was pointed out by Bose and Einstein in 
1925. On the level of mathematical physics, the understanding of the phase transition 
started with the well known paper of Araki and Woods [1], where the generating func- 
tional of the cyclic representations of the canonical commutation relations corresponding 
to the equilibrium states of the free Bose gas are computed for periodic boundary con- 
ditions. Lewis and Pule [7,8,11] computed the grand canonical equilibrium states for a 
set of boundary conditions including the Dirichlet and Neumann boundary conditions but 
not the attractive boundary conditions. They are using the Kac method. An important 
consequence of their result is the explicit computation of a non-trivial Kac density showing 
non-equivalence of the canonical and grand-canonical ensembles in the condensate region. 
The next result is found in [3], where the same conclusion was obtained for generalized 
condensations in some models of imperfect gases with diagonal interactions. 

In the present paper we complete this computation of the equilibrium states for the 
free Bose gas with attractive boundary conditions. About the relevance of this type of 
boundary conditions, see e.g. [4, 10]. This model has a particular type of condensation 
namely condensation in quantum states corresponding to isolated points in the spectrum. 
It is well known [10], [5] that in this case the condensate is situated at the "boundary" 
and not uniformly spread out everywhere in space. We give a precise formulation of the 
generating functional in the frame of the theory of generating functionals on the CCR in 
order to catch up the condensate. Finally we derive also that there is non- equivalence of 
ensembles, something which was unclear until now because of the fact that the quantum 
fluctuations show a pattern [6] completely different from the free Bose gas with Dirichlet 
or Neumann boundary conditions. The intuition behind this fact is related to a wondering 
peculiarity of the free Bose gas with attractive boundary conditions [10], [5]. If one takes 
the thermodynamic limit using a sequence of growing convex domains with the point of 
homothety at the origin of the coordinates, then locally the condensation density is always 
equal to zero. As a byproduct, our computations imply that the condensate is spatially 
situated in a region logarithmically close to the boundary of these increasing domains. In 
the present paper we take different positions of the homothety point for cubic containers 
to show that the density of this logarithmic stratum of condensate inherits also a spacial 
anisotropy due to the choice of cubic containers. Finally remark that for the rotating 
bucket case [11], one has also the effect of the condensate being increased at the boundary. 
But this is an effect of large angular momentum and not of the boundary conditions as in 
our case. 
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2 CCR- Representations and the generating functional 



For details about the CCR algebra, we refer to [2]. 

Let h be a complex pre-Hilbert space with inner product (•,•). A representation of the 
CCR over h on a Hilbert space H is a map / h- > W(f) of /i into the group U(H) of unitary 
operators on a Hilbert space satisfying the Weyl relations: 

W{h)W{f 2 ) = exp j-i Im^, / 2 ) } W{h + / 2 ) (2.1) 

such that for each f Eh the map A i— > W(A/) of K. into U(7i) is strongly continuous. By 
Stone's theorem, this continuity condition implies the existence of self-adjoint operators 
$(/) such that 

= exp{id>(/)} (2.2) 

These $(/) are called field operators. The map / i— > $(/) is linear over R, but not linear 
over C. Using the $(/) we can now define the creation and annihilation operators a*(f) 
and a(f) for f E h by: 

a*(/)=2- 1 /2 { $( / )_ i $( i/)} (2.3) 

a(/) = 2- 1 / 2 {$(/) + i$(i/)} (2.4) 
A state on the CCR-algebra is a linear functional u> : h — > C with the properties: 

cj(A*A) > 0, u;(l) = 1, for A linear combinations of the f E h 

A representation (W, 7i, Q) is called a cyclic representation if Q is a cyclic vector. A vector 
fl is cyclic if the set {W(f)Q} is dense in 7i. To each cyclic representation (W, TC, Q) 
of the CCR corresponds a generating functional E : h — > C given by: 

E(/)=c(H/(/)) = (fi,^(/)fi) (2.5) 

Proposition 2.1. ^4 functional E : /i — > C £/ie generating functional of a cyclic repre- 
sentation of the CCR if and only if it satisfies the following conditions: 

(i) E(0) = I, 

(ii) V/ E h : A i— > E(A/) continuous , 

(iii) V /mite sets o/ complex numbers c±, . . . , c„ and elements fi, ■ ■ ■ , f n E h : 
E,E,E(/,-/,)ei Im ^^c- Cj >0. 
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3 Kac-density and equivalence of ensembles 
3.1 Concrete setup 

Let A£ = [— L/2,L/2] !/ be a bounded region in W with volume V = U . We put hi = 
C 2 (A U L ) for the Hilbert space of the wave-functions in A U L with the scalar product (/, g) h : = 
J* A „ dx u f(x)g(x). Then K V L C A£, and C whenever L < L' via natural imbedding. 

Let be the self-adjoint extension of the operator — (with domain dom{— A^) = 
C^°(A£)) determined by the boundary conditions d n (p + <t0 = on <9A£. Here d n is the 
directional derivative in the direction of the outward normal n to <9A£. If the parameter 
a < 0, we say that the boundary <9A£ is attractive. 

First we have to solve the one-dimensional one-body eigenvalue problem onA L = [-L/2, L/2]: 



Due to these attractive boundary conditions, there are two negative eigenvalues tending 
to the same limit —a 2 (when L — > oo) and an infinite number of positive eigenvalues (for 



if L 4>) (x) = \4>(x) 



with boundary conditions (a < 0): 




L\a\ > 2): 



e L (0)< ei (l)<0<e L (2)< eL (3)<..., 
C£ (0) = -<7 2 -O(e- L W), 
c L (l) = - < 7 2 + 0(e- L W), 




(3.1) 



The corresponding eigenfunctions {4>k}kei+ form a basis in h L and are given by 




for k even, 



for k odd. 
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The eigenvalues and the wave functions of the corresponding multi-dimensional case have 
the form: 

V 

E L (k) = J>*)> 

V 

v£(x) = hum, 
i=i 

where k = {hY i=1 G 7L\ and x = {xi} u i=1 G k v L . 



3.2 Kac density 

The Kac density relates expectation values of observables in the canonical ensemble and 
those in the grand canonical ensemble. The canonical equilibrium state for a free Bose gas 
in a cube K V L of volume V — If with total particle density p and inverse temperature f3 is 
given by 

T V „) AW e-^ (n) 

"SUA) = — ^ -sag-, where n = [Vp], M^') C ?#\ (3.2) 

Tr („) e ^ 

and T^'^ is the n— particle free Bose gas Hamiltonian in the cube K V L with boundary 
conditions defined by a. Now we consider the grand canonical equilibrium state at chemical 
potential p and inverse temperature (3: 



W L,/3,aA A ) = — -— r o/™ —r M ' rfom(A) C J^L,B- (3.3) 



Tr^ |B Aexp{- j 3(Tf -/xiV £ )} 
Tr^ s exp{-/5(T£-/iiV L )} 

Here T£ = Xlkez^ -^i(^) a *(V ; k) a (' ? /'k) i s the f ree Bose gas Hamiltonian and N L = J2kez^ ^L,k 
Skez" a *(V ; k) a (' ? /'k); i s the particle number operator in Tl,b, the boson Fock space over 

Tl,b = Fb(C\V l )) 

oo 

= ©Wfi (3-4) 

n=0 

with the symmetrized n— particle Hilbert space appropriate for bosons and Tt^ B = C. 
Notice that in the thermodynamic limit L — > oo the canonical ensemble state u/g a ™(-) may 
not coincide with the equilibrium state of the grand canonical ensemble state Vp^m p )( m ) 
for the corresponding particle density p. Here ~p(/3,p) = lim^oo ~p L ((3, p) and ~Pl(P,p) is a 
solution of the grand canonical particle density equation (see also (|3.22|) ) 

P = "LM M (N L /V). (3.5) 
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By virtue of ()3.4|) the states ()3.2|) and ()3.3|) are related by 



XA) = I K'fijtf) u^(A m]) l (3-6) 

JR + 

where = A\Tl£ b is a restriction of the operator A on the subspace 7v£ B and 

ES ex P ( n ^) ^^1% exp ( "TO ) 
= Tr^expj-^-^)} • (3 - 7) 

For a given grand canonical density (|3.5|) . the measure f)3.7j) takes the form : 

JiW*) : = K9 Lh L (M^ = ^K L ^p), (3.8) 

and the limit Kp(x; p) = lim^oo Kl^(x; p) is known as the Kac density, see e.g. [8]. If the 
Kac density happens to be a 5-function with support at p , then clearly one has [strong) 
equivalence of ensembles: 

"t^M)=bftW> ( 3 - 9 ) 
Otherwise there is only weak equivalence of ensembles, see [3]. 

The limit p c (P) '■= lim^-^ lim^oo w 9 £ C p (Nl/V) is the critical density for the free 
Bose gas in a box with attractive boundary conditions. We shall show that in the model the 
canonical and the grand canonical ensembles are not equivalent in the presence of the Bose 
condensate, i.e. for p > p c (P), or for (3 > /3 c (p), where p c (/3 c (p)) = p. The non-equivalence 
of ensembles in the case of the free Bose gas with attractive boundaries is not the same 
phenomenon as in the case of the one with, for example Dirichlet, a = oo, or Neumann, 
cr = 0, boundary conditions. In the case of the attractive boundary conditions (er < 0), 
the condensation phenomenon is a surface effect (not a bulk effect as in the free Bose gas 
with a = oo or o = 0): the condensate is located near the walls, see Section 4.2. 

To determine the Kac density, we have to calculate (see (|3.6j) ) 



d^K L ^;p)uTUW(f)) (3.10) 



for any test function / G Cq°(W), the C°°-functions on W with compact support. There- 
fore we first must calculate the limit of the expectation value of the exponential function: 

with {p L (fl, p)}l solutions of the density equation (|3.5|) . 



6 



This is possible because the states u 9 £ C p , where \i < —ucr 2 , are quasi-free states, and 
these are easily obtained by using the truncated junctionals 0JL,p,p{- ■ -)t, see e.g. [2]. The 
functionals are defined by the recursion relations: 

<A^i ■■■A n ) = E II <P,Mm> • • • » A W\)h (3-12) 

t£P„ JSt 

for all Ai (i = 1, 2) creation or annihilation operators and neR The sum r G P n is over 
all partitions r of a set of n elements into ordered subsets J = {j(l), . . . , j(\ J |)} G r. 
One can verify that the truncated functionals associated to the equilibrium states u 9 L '°* 
satisfy 

<J>*(/))t = <£>»(/)) = 0, 

<A>*(A), «*(/*))* = ^(a(/i),a(/ 2 )) r = 0, 

<A>*(A).«(/a))T = f/2, /i) , (3-13) 

V e ^ f L _1 / ^ 

with f,fi,f 2 ,... G the space of testf unctions with support in A£, a" = {a or a*} and 
is the self-adjoint extension of the Laplacian — corresponding to attractive boundary 
conditions o < on dh v L . Then the non-trivial two-point functions ()3.13|) are explicitly 
given by 

^5>*(/i),«(/ a ))r = «£5>*(/i)a(/ a )) 

= EAOOAW ^w-m)-! ' (3 ' 14) 

where the transformation /(x) t— > /(k) of / G Qj >c (IR ly ), is now defined by 

/(k) := (^, = / rfx ^JM/(x), (3.15) 
JaI 

the Fourier transforms for the basis of t a L (see section l3~T|) . 
Now 

oo . n 

= expV^ ij/3iP ($(/), $(/),..., $(/)) T , (3.16) 
^— f n! > v ' 

n=l , . 

n times 

where the ($(/), <&(/), . . . , $(/))t are the n— point truncated field correlation func- 
tions. Because of the fact that w£'^ iAt is a quasi-free state, only the two-point truncated 
correlation function is non-vanishing, yielding: 

rtsJWti)) = exp (^J* W), *(/))r) (3.17) 

7 



By virtue of ()2.3|) and ()2.4|) it can be rewritten in terms of the creation and annihilation 

operators a*(f) and a(f): 

$ (/))^ = \<J> /K + <S>*(/M/)) (3-i8) 

so that the explicit form of the generating functional ()3.17|) becomes: 
^(W(f)) = r ^^(f(/)) 

fir- L— >oo ^ 

= e^^-i(/,/) fe -i2™w£5>*(/)«(/))) ( 3 - 19 ) 

A last remark about the thermodynamic limit. Notice that the grand-canonical ensemble 
for the free Bose gas exists only for p < inf spec(t u L ). Therefore, the solution of equation 
(13. 5|) verifies the inequality p L (j3,p) < —ua 2 . Since the critical density 

lim lim <^ = p c {(3) (3.20) 

for the free Bose gas with attractive boundary conditions a < is finite for all dimensions 
greater then, or equal to one, [10], [5], Bose-Einstein condensation occurs for p > p c ((3)'- 

p {(5) := p - p c {(3) = Km r^- L(M (^j > 0, (3.21) 

where Nl,o is the number-operator on J~l,b of the zero mode k = 0. The factor 2 V is due to 
the asymptotic degeneracy of the inf spec(tj) = — vo 1 + 0(e _L ' <7 ') for L — > oo , see Section 
3.1. Notice that (|3.21|) implies that the solution of (j3.5|) for p > p c (P) has the asymptotics: 

MM = - v S- w -^ m - + 0{ L-»). (3.22) 

We use this result in the computations of the thermodynamic limit of the generating 
functional below. 

We conclude this section by the following statement about the explicit form of the Kac 
density for the thermodynamic limit of the free Bose gas in the cubic box with attractive 
boundary conditions. 

Theorem 3.1. For the free Bose gas T[ with attractive boundary conditions a < the 
limiting Kac density has the form : 

Kpfop) (3.23) 
8(£-p), forp<p c ((3), 



2"fl(g-p c Q3)) 
(2"-l)!(p- Pe C9)) 



P - Pc{P) 



2"-l 



cxp 



rwr- >, for p> p c (P). 

P ~ PcKP) J 



Here 6{z < 0) = and 6{z > 0) = 1. 
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Proof: By the identity (|3.10p . the Kac density Kp(£;p) is related to the thermodynamic 
limit of the characteristic function of the particle density N^/V for t G K 1 : 

lim 

L— >oo 



/ de^(e;p)ex P (zte), (3.24) 



To calculate the limit in the left-hand side of (|3.24|) . we use that the state u 9 L C p-,p ^(-) is 
quasi-free. Then 

<^(Ap)(e*P(^V)) 

TT f l-exp[-/?(E L (k)-/I L (Ap))] 1 
k ^ll-exp[-/3(^(k)-/7 r (/3,p)-zt//3^)]|- l - J 

Since the 2 U lowest energy-levels, i.e. the levels for which k G K<2^ = {k G : hi = 
0,1; z = 1, . . . , v} are exponentially degenerated when L — > oo : ^(k G K< 2 ^) = — z/o" 2 + 
0(e" L|ff| ), by virtue of (IT2^|) and (ET231) we get that 

^L% L v>, P )(?xp(i tN L/v)) 

exp(ztp), for p < p c (/3), 

[1 _ zt2 ^(p - Pc (/?))r 2 ' exp(z*p c (/3)), for p > p e (J3). 1 ' J 

Therefore, by (|H.24jl . the Kac density (|H.2Hjl is the Fourier transformation of the right-hand 
side of (ET2EI) . □ 



4 The generating functional 

4.1 Condensate and generating functional 

We are interested in the thermodynamic limit of the generating functional Up C p {W(f)) = 
liniL-ooa;^^ ^ } (W(/)) for any / G CE°(R"). To this end we choose the box k u L with 
L large enough such that the A/ := supp(f) is contained in A£. We consider here the 
generating functional Up„(W(f)) for / an element in Cq 3 (M u ). 

Theorem 4.1. The generating functional uip C p (W(f)) on C^°(R U ) is given by: 

^ P (W(f)) = ex V (~(fJ)y w (~(f,g a ((3,P)f)), (4-1) 



9 



with operator g a (/3 , p) on £ 2 (R U ) defined by 

(a,OM/)(x) = I dy G CT (Ap)(l|x-y||)/(y), 

e -npp(J3,p) 



G a (P,p)(r) 



(4ir(3)-»/ 2 J2 e ~ r2/ ^ 



n=l 



n 



v/2 



(4.2) 



where p(P,p) < —vo 2 for p < p c (P) and p(fl,p) = —vo 2 for p > p c ((3)are limiting solutions 
of the grand canonical density equation i3.5}) . 

Proof: In order to determine the generating functional Wp°„[W(f)), we have to compute 
limL_ > oo<^£'^ i ( / gp)(a*(/)a(/)), see (13.19)) . Since the attractive boundary conditions o < 
create a gap in the spectrum spect(t1), and respectively in spect(T[), the calculations 
need a separation of the negative eigenvalues from the positive part of the spectrum. 

We consider first the one-dimensional case, when there are only two negative eigenvalues 
tending to —a 2 for L — > oo, see Section 13.11 By virtue of ()3.14j) one gets for a given 
/ E C^iR 1 ) that 



lim 

L— >oo 



< c ^ L (/3>*(/M/)) 



lim l/>)P 

L— >oo ' * 



1 



1 



j!S,(l/( )l' e««W-^^)-l +l/(1)|1 

OO oo 



3 /3(e i (l)-/I ij (/3,p)) _l 



+ ^ e n^(/3,p)^ e -n/3e iW | /( ^ | (43) 
n=l fc=2 

As mentioned before, we choose large enough such that supp(f) = Af is contained in 
A^. Then one estimates that |/(0)| 2 has an asymptotics of the order of 0(e~ i ' cr ') for large 
L since 

2 



l/(o)P 



2 / + sinh(LH) y 1 /■ rfa;/(x)cosh( _ Hx 

L \ L\<7\ J JA f 

/ dxf(x) cosh(— \a\x 

J At 



L 

4|a|e- L H 



+ o e" 



(4.4) 



The integral in the last expression is independent of L, because supp (f) is finite and inside 
the box A^. Similarly one gets for L — > oo that 

2 



|/(l)| 2 = 4|a|e-^l 



dxf(x) sinh(— \o~\x) 



+ o e 



-L|of 



(4.5) 
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Consider now the coefficients of |/l(0)| 2 and of |/l(1)| 2 in (jOJ). If p < p c {ft), then Jt(/3, p) < 
-cr 2 , i.e. ei,(0) - ~p, L (fi,p) > for large L. Therefore, by virtue of (fO)l . (|Q)l and (|Oj) . 
both of those terms are of the order 0(e _L ' CT ') for large L. If p > p c ((3), then JI L (f3,p) = 
—a 2 + O^L" 1 ), and one gets for large L: 



3 /3(e L (0)-/I L (/3,p)) _x U V ^ e 9(e i (l)-7* x ,08 I p)) _l 

^^o(/3)l{|/(0)| 2 + |/(1)| 2 }, (4.6) 

where po((3) = p — p c (/3) is the condensate density. Therefore, again by virtue of ()4.4j) and 
(J4.5|) . these terms vanish in the limit L — > oo. 

Consider now the last term in the limit (|4.3|) . By virtue of (|3.15|) for u = 1 (see Section 
3.1) we can represent the sum over > 2 in the following explicit form: 

2 ^ f i + MVW)L) \ 1 j^r )heUk) e -^ (fc) + 

1 { k>2: even} \ V^(k)L J 

L 



{k>2: odd} 

where s = n(3 and 



2 E / ^(^^ N- 1 ^-^ ^ 



(4.7) 



FLeven(k) := / dx cos(^fTJJz)x) f (x) , F Lodd (k) : = / rfx sin( A/e L ( / (x) . (4.8) 

./A/ ^A/ 

Since the spectrum {ei(k)} k>2 verifies the conditions ()3.1|) and / G (^(M 1 ), the first 
and the second series of terms in (|4.7jl are Darboux-Riemann sums for the corresponding 
integrals: 




1 + sin (vgg L ) ^| F Leven (k) F Leven (k) e-*« 



{k>2: even} 



V^ik)L j 



+ 2 / sin(y^(fc)L) 



F Lodd (k)F Lodd (k) e— < fc > 



- J dfcHe(e»-,/) fe Re(e 4fc ,/) / 

i roc 

+ - dk lm(e^-,f) hL Im(e lk -j) hL e- sk2 . (4.9) 

J 



The last expression of (|4.9J1 yields 
1 



2tt / v ^ ' J "* v ^ ,JJh ^ 



dk(e^J) hL {e ik ,f) t v 



{Aits)- 1 ' 2 I dx I dyf(x)f(y) exp <| L . (4.10) 



4s 
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Finally, taking into account (|4.4| ) -([4.6 |) . (|4.1U|) . and the fact that p(f3,p) < —o 2 < 0, we 
get for the limit ()4.3|) in the one-dimensional case: 

"£>*(/)«(/)) = (f,g^=i(P,p)f) hL , (4.11) 
where g a ,u=i(P,p) is the integral operator on £ 2 (M 1 ) defined by 

(g*,u=i(P,p)f)(x) = / dy G a>u=1 (P,p)(\x -y\)f(y) , 

e -npp(J3,p) 



G^ =1 (/3,p)(r) = (4 7 r/3)- 1 / 2 J]e 



-r 2 /Anf3 



n l/2 
n=l 



Using the results for the one-dimensional case, one computes the two-point correlation 
function u 9 ^ p (a*(f)a(f)) in the //-dimensional case. Since the first 2 U wave functions 
(x)} kgZJ/ have the same exponential behaviour as in the one-dimensional case and 
since inf L special) = —ua 2 , see Section 3.1, we get: 



2™ </3,P(/3,p)( G *(/) a (/)) 

\fW\ 2 e f3(E L (V-Tl L (M)-l 

OO p p 



71=1 



4n/5 



which implies (j4.2|) . By virtue of f|3.19|) this finishes the proof of (J4.1j) and of the theorem 
for any particle density p. □ 



Theorem 14. II tells us that the condensate is not traceable by considering only strictly local 
observables. The characteristic functional on the CCR-C*-algebra of quasi-local observ- 
ables coincides with the one without condensate. The reason for this is that the condensate 
is not homogeneous but located in the vicinity of the container boundary. 

In order to catch up the presence of the condensate or to get a complete picture of 
the system, one has to extend the algebra of observables to the weak closure of the CCR- 
C*-algebra with respect to the limit Gibbs states. In the next paragraph we compute the 
limit functional on the relevant non-localised observables, and obtain a complete picture 
yielding the existence of sufficiently many fields in the representation of any w*-limit point 
of Gibbs states as L tends to infinity. In fact our strategy will be to make a relevant choice 
of the homothety point for the thermodynamic limit of convex containers, in order to catch 
up the condensate. 

Above and below we considered only the easy shape container limit, namely cubic boxes. 
Because of the particular inhomogeneous spreading of the condensate in the neighbourhood 
of the box boundary, it is clear that this thermodynamic limit treatment can be very much 
shape dependent. In this paper we do not enter into the details of this specific problem. 
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4.2 Condensate localization 



Remark 4.2. It sounds curious that in spite of the non-zero condensate density for p > 
Pc{/3) ; f^S.21\) . there is no trace of it in the generating functional \4-l\) . This is in contrast 
to the Kac density hS. 2ty . which explicitly depends on the condensate density p — p c (f3). To 
understand this difference one has to take into account that fl^.i| ) is localized on the support 
of the function f e C£°(M. V ) whereas the Kac density is a global function, depending on the 
condensate even if it is localized at "infinity", sticked to the attractive boundaries. 

In order to make this statement rigorous we start first with the one-dimensional case. 
Let the function / e C^M 1 ) be such that supp(f) = (-6,5) C (-L/2,L/2) and 5 < 
(lnL)/2|a|. Consider its shift over a distance 7l(ct) := L/2 — (2|cr|)^ 1 InL: 



KtM = :=f{x- [L/2 - (21a!)- 1 InL]) 



(4.12) 



Then f TjL(a) GC -(-L/2, L/2). 



To get the generating functional we compute now the limit of the corresponding two- 
point function (|3.19jl : 



lim u 



L— >oo 



g. c. 

L,f3,p L (p,p) 



ax 



7LW/) a ( T 7iW/)) 



t^^^^)^ 2 e f3(e L (k)~JT L (M) 
k=0 



lim \f T . 



7lO) 



(o)p 



3 0(e L (O)--p L (l3,p)) _l 



+ \fi 



7l( ct ) 



zt+EiW*>i ! 



1 



k=2 



3 (3(e L (k)--p L ((3,p)) _l 



.(4.13) 



Remark 4.3. Notice that in contrast to \4-ty , the shift \4-12]j corresponds simply to the 
choice of a new point of homothety for the thermodynamic limit \4.13\) . In \4-3j} , the point 
of homothety coincides with the origin of coordinates x = 0, whereas in \4-lty this point is 
L/2 - (2\a\Y l \nL. 

Now, and in contrast to (|4.4|) . |/ t ^ (ct) (0)| 2 goes like L" 1 for large L. Indeed, 



I/t 7l(ct) 



1L (?)+& 



\o-\L 



-i 



dxf{x)e^ x 



+ o{L 



Remark that for p > p c {/3) the first term in ()4.13|) remains now finite in the limit L — > oo. 
Taking into account ()3.21|) and ()3.22|) one gets: 



lim |/ T , 



7l(o") 



Po{P,p), 



*0(eL(P)--p L (p,p)) _1 



(7 



5 2 
\a\x 



dxf(x) &' 



(4.14) 
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The same reasoning for the second term in formula (j4.13j) gives a similar result: 

1 Po(/?,p) 



\CJ\ 



5 2 
\er\x 



dxf(x) &' 



(4.15) 



By the same computations as used in the proof of Theorem 14. 1[ the third term in ()4.13|) 
yields for p > p c (ft)'- 



i 



k=2 

oo 



e P(e L (k)-]I L (l3,p)) _l 



Ve- n/3ff2 (4vrn/3)- 1 / 2 f dx I 



dyf{x)f(y)exp 



\x - y\ 

An(3 



(4.16) 



Hence the two-point function for the one-dimensional problem becomes: 



i^o W S,7I L (/3,P) ( fl *( T 7L (<x)/)«(^l(-)/)) 

2 

dxf{x)e^ x 



Po(P,p)W\ 



+ (f,ga,u=l(P, P )f), 



see (|4.11|) for the definition of the operator g a ,u=i(P, p)- 

It is evident that one gets the same result for the shift of supp(f) 
distance -Jl(ct) = -L/2 + (2\a\y 1 InL, i.e.: 



(4.17) 
-S, S) over a 



= Po(/3,p)M 



a*{r± 1L (a)f)a{r± 1L {a)f)) 



dxf(x) e 



±|cr|x 



+ (/,V=i(Ap)/) 



where 



(T± lL (<r)f)(x) := / (x=F [L/2 - (21a!)- 1 InL]) . 



(4.18) 



(4.19) 



Therefore, taking the thermodynamic limit L — > oo at one of the homothety points ±7^(17), 
we get that the generating functional depends on the Bose-condensate density for p > p c (/3): 



= exp (-\(f,f)\ exp Lhj% =l {f) - \{f,g„,v=xf)} , (4.20) 
where, by virtue of ()4.18|) . one has 



Ct v=1 (f) = Po{P,p)\a\ 



dxf(x) e 



± I crl a: 



(4.21) 
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Remark 4.4. Notice that this result is due to a fine (logarithmic) tuning of the position 
of the homothety points ±72,(0"). Indeed, take ±72(00), for < a < 1, i.e. the homothety 
points are more distant from the boundary ±L/2. Taking into account the explicit form 
of the eig en functions for k — 0, 1 one finds that now \ f T±n/ (aCT) (A; = 0, 1)| 2 goes for large 
L like L~ l l a . This implies that both limits \4-14\) and \4-15\ l, and hence \4.21^ , vanish. 
So, the generating functional \4-2(f^ has the same form as for thermodynamic limit with 
the homothety point at the origin. In contrast to that, the choice 1 < a means that the 
homothety points are closer to the boundaries ±L/2. Then |/r ±7i(a£r )(^ — 0, 1)| 2 goes slower 
then Lr 1 . This implies that both limits \4-14\j and and hence \4-21\ j, becomes 

infinite. So, for p > p c (/3) the generating functional \4-%Oj) zero, whereas for p < p c (/3) 
it is nontrivial with C av=x {f) = 0. 

To interpret these results, consider the local particle density: 

fcez^ =1 

Here a(x) is the Bose-field operator such that a(f) = L u =i dxf(x)a(x) for / e C^°(A^ =1 ) 
and N(x) = a*(x)a(x) is the local number operator, cf. ()3.14|) . Then by (j3.5|) and f)3.14|) . 
the global density is: 

<ln L (M (x) = I j Arl dxuJ W L tM ( a *( x W x )) ■ ( 423 ) 

Consider the thermodynamic limit of the local particle density at the origin of the coor- 
dinates x = 0. Taking into account the explicit form of the eigenfunctions, one gets that 



P(P, P ; * = 0) := hm <^ (Ap) (a*(x = 0)a(x = 0)) = \ dk ^_^ p)) - (4.24) 
for p < p c ((3), and 

If 1 

p(P,p; x = 0)= lim <^ ()M (a*(x = 0)a(x = 0)) = - / dk /3(fc2+fj2) _ = p c (P) 

(4.25) 

for p > Pc(ft) by ()3.22j) . By inspection of ()4.24j) and ()4.25|) based on the explicit formulae 
for the eigenfunctions one readily gets that 

P (f3,p;x)=p(f3,p;x = 0) (4.26) 

for any 1 in a bounded domain D, containing the origin of the coordinates x = 0. In 
particular we get that the limiting local density for x G D corresponding to the first two 
modes {k = 0, 1) is 

p ((3,p; x) := lim V -Mfeff^r = 0. (4.27) 

fc=0,l 
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On the other hand, the global Bose-Einstein condensation density (j3.21j) is also related 
exactly to these two modes: 

Po (P, p) = lim \ £ e ^ L[k) J, LiM) ZI = P ~ P*W > 0. ( 4 - 28 ) 

fe=0,l 

which is not present in ()4.25|) . 

Consider now the local density of the Bose-Einstein condensation ()4.27|) at the homo- 
thety points ±7i(cr). Then taking into account the explicit form of the eigenfunctions 
(j>k=o,i{ x ) an d (|3.22j) . we get that, in contrast to (|4.27j) . the local condensate density is 

lim V l fl:-l L{ ^ -Po(P,p)\a\. (4.29) 

fc=0,l 

The same arguments as above show that this condensate local density varies from zero 
to infinity when the parameter a in the homothety point positions ±7i(acr) varies in the 
same interval. 

Remark 4.5. These observations can be interpreted as follows : the Bose-Einstein con- 
densate for attractive boundary conditions a < is localized in a logarithmically narrow 
domain in the vicinity of the boundary. In other words this kind of condensation is a 
surface phenomenon. At the same time globally it is very "visible", since the Kac density 
indicates a non- equivalence of ensembles in the presence of the condensate, see Theorem 
3.1. 

For the generalization to the i/-dimensional case, we start with the corresponding local 
condensate density: 

p (/3,p;x):= hm eflj££ol)) -1 ' (430) 

ke{Z^: fe a =0,l ; a=l,...,v} 

Let x belong to a bounded domain D u , containing the origin of the coordinates x = 0. 
Then using the explicit expressions for the eigenfunctions ^k( x ); see Section 3.1, and by 
the same arguments as above for v = 1, we obtain that the limit ()4.30|) is zero for all 
densities p > 0. 

The product structure : ^k( x ) = T[i=i ( l > k ( X/ di hnplies that this conclusion does not 
change if we consider instead of x £ D", the condensate density in the vicinity of the points 
corresponding to the shifts where at least one among the v arguments remains unshifted. 
On the other hand, this structure and the asymptotics of (p^ =0 x(x) for \x\ — > oo yields also 
that for any k £ \7h\ : k a — 0, 1 ; a — 1, . . . , v} one gets 

(il r ±T*(-* ( x = °) = | ( rr/ 2 L-K«r 1 +« 2 - 1 +-+^ 1 ) + o (VH^+^+.-.w)) (431) 
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as L — > oo. Then, by virtue of (|3.22|) . the limit for the local condensate density becomes 
non-trivial: 

1- I (n«=l r ±7£(a-aO-)V'k) ( X = 0)| I if / ^ n qo\ 

&, 2s e p(E Li ^ L (M) _i = M (p-Pc(/3))>0, (4.32) 

ke{z^:Jfc a =0,l;a=l,...,i/} 

if and only if 

af 1 + a 2 x + . . . + a" 1 = z/ . (4.33) 

This means that the condensate (up to logarithmic deviations) is localized essentially in 
the corners of the hypercube A£, where L — > oo. We proved the following statement: 

Theorem 4.6. Let li. be in a bounded domain D u , containing the origin of the coordinates 
x = ; then the thermodynamic limit of the local particle density is 

p(P, p; x) := Km £ ^ ( J^j, p)) = ^ / ^(k^L)) - , (4.34) 
kez^ ' /K + 

where Jl({3,p) < —ua 2 for p < p c (/3) and~p((3,p) = —ua 2 for p > p c {(3). Thus p(/3,p; x) = 
p c (P) for p > Pc(P), i-e. the local condensate density p (j3,p; x) = for any p > 0. 
Whereas at the homothety points corresponding to the shifts Yla=i T ±-iL(a a <r) with parameters 
satisfying \Jf..3$) , the local condensate density H4-3£\ ) is nontrivial. Moreover, besides being 
inhomogeneous it is also anisotropic and essentially localized in the directions of the corners 
of the hypercube A^^. Varying the parameters {a a } v a=1 in the range (0, +oo) one finds 
this local condensate density varying from zero to infinity. 

Now we extend Theorem 14.11 on the basis of our discussion above of the condensate lo- 
calisation and Theorem 14.61 Similar to the one-dimensional case, see Remark 14.41 our 
relevant localized observable in the i/-dimensional case will be a function / 6 C£°(W) such 
that supp(f) = (—5i,5i) x (—62,62) x ■■■ x (—6 V ,6 U ) C K U L and with 6 = max i=lj jjy 5i 
such that 6 < (lnL)/2|cr|. Consider in each coordinate the shift over a distance 7z,(er) = 
L/2- (2|a|)- 1 lnL: 

f[ r jUa) A (x) =/(*!- (L/2 - (21a!)- 1 In L), . . . , x v - (L/2 - (2M)- 1 InL)) , 

\a=l J 

thennLi^(<r)/eCo°(A£). 

To get the generating functional in the i/-dimensional case, we compute the limit of the 
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corresponding two-point correlation function: 



Yl T 7L^)f a Yl T jL^)f 



\a=l 



\a=l 



lim > 



lim 

L— >oo 



E 

ke{Z^:fc a =0,l;a=l,...,^} 



2 



a=l 



a= 

E 



1 



d P(E L (k)--p L (f3, P )) _l 



n=l 



a=l 



(4.35) 



k€Z^\{Z£: fc a =0,l; a=l, v} 

This thermodynamic limit depends on the homothety point corresponding to the shifts 



n«=i T ±-y L {a a <r) with parameters a a = 1. Notice that the factor Yla=i /^mOO 
order 0(L~ U ) for large L: 



is of the 



n ^ L ( CT )( k ) 



7 L (cr)+(5i 

7 L (o-)-5 
li/ r — u 



dx v n^w^w^M 

Ti(o-)-^ q,=i 
2 



supp(f) 



dx/(x)Hel' 



+ o(L-"), 



(4.36) 



and for any k G {Z^ : fc a = 0, l;a = 1, Hence, by the same reasoning, which 

implies (|4.32j) . the first 2 U terms in (|4.35j) give: 



lim 

L— >oo 



E 



ke{Z^:fc Q =0,l;a=l,---^} 



II(^>)( k ) 



1 



po(/3,p)M 



dx/(x)ne |CT|si 



(4.37) 



For the last term in ()4.35|) . we perform the computations as in Theorem 14. 11 yielding: 

2 



lim 

L— >oo 



E 



a=l 



1 



3/8(£zi(k)-TiiGM) -1 



keZ^\{Z^:fc a =0,l; a=l, *4 
00 

/^ (47rn/3 )-/2 / dx / rfy/ (x)7(^)exp 



n=l 



l x ~y| 



4nf3 



(4.38) 
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So, taking the thermodynamic limit L — > oo at one of the homothety points {±7l(<t)}^ =1 , 
we get now the generating functional for p > p c (/3): 



cxp 




■ S (/,/))exp 



■7;C£(P,p){f)--(J,ga(P,p)f) 



with 



C±(P, P )(f) = p Q (P,p)\<r\ v 



dx/(x) JJe ±|ff|a;D 



a=l 



(4.39) 



(4.40) 



Remark 4.7. Again this result is due to a fine (logarithmic) tuning of the position of the 
homothety points {±7^(cr)}^ =1 in the corner directions of the hypercube K V L . Indeed, take 



as in Theorem 



the shifts Y[ a=l T± lL{aaa) , with ]T 



a=l 



> v, i.e. the homothety 



points are more distant from the corners of the hypercube. Taking into account the explicit 
form of the eigenf unctions for k 6 {Z 1 ^ : k a — 0, 1; a — l,...,z/}, one finds now that 
I Ua=i fr JL(aa ^)\ 2 with k 6 {Z* : k a = 0, 1; a = 1, . . . , u] goes like L-KW+.-W) for 
large L. This implies that the limits of the first 2 U terms \4-3 r i\) , and hence \4-4$ > vanish. 
So, the generating functional has the same form as for the thermodynamic limit with 

the homothety point at the origin x = 0. In contrast to that, the choice < Yla=i (°a) _1 < 



v means that the homothety points are too close to the corners of the hypercube K U L . Then 



Yla= fr±-y T ( a „a)0 i -)\ 2 k e {^+ • — 0,1; a = l,...,z/} goes to zero slower then 



L v . This implies that the limit becomes infinite. So, for p > p c {/3) the generating 

functional J^.ff£| ) is zero, whereas for p < p c {(3) it is nontrivial with C^(f) = 0. 

Therefore, we proved the following theorem: 

Theorem 4.8. The generating functional lim^. 
is given by 



^LATX L ( M {U: =1 W(r ± , M f)) onC^il 



\[W{r ±<aaa) f) 
1 



va=l 



exp ( ~^(f,f) ) exp 



cttf,p){f)-M9Mp)f) 



(4.41) 



with 



c£GM(/) 

(g a ((3,p)f)(x) 



Po(P,p)W\ 



dx/(x)JJe ±k|! 



i=i 



dyG a ([3,p)(\\x-y\\)f(y), 

-n/3va 2 



-r 2 /4nf3 



n=l 



n 



v/2 
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Here JZ(/3,p) < —vo 2 for p < p c (/3) and ~p(j3,p) = —va 2 for p > p c (P) are the limiting 
solutions of the grand canonical density equation \3. <5j) . po((3,p) = for p < p c (/3) whereas 
Po{fl,p) = P — Pc{0) f or P > Pc(P) and the function Xv{ a ) — 0, 1, +oo, respectively for 
v < ELi (a-a)' 1 , v = E«=i M' 1 , and v > J2a=i W" 1 - □ 

5 Concluding remarks 

The main results of our analysis for the free Bose gas with attractive boundary conditions 
are contained in the Theorems 13. 1| 14.11 and 14.81 

In Theorem 13.11 we obtain a Kac density function showing non-equivalence of the 
canonical and the grand canonical ensemble in the presence of the condensate even if the 
condensate density is locally zero. 

We learn from Theorem l4. ll that the condensation is not visible in the expectation values 
of strictly localized observables because the Bose condensate is situated near the boundary 
of an "infinite container". Nevertheless one should observe the effects of condensation in 
the equilibrium states, i.e., in the generating functional. 

Theorem 14.81 yields the answer. We make precise for which type of observables the 
equilibrium states show their dependence on the condensate. This completes the rigor- 
ous analysis of the problem of (non-)equivalence of ensembles for the free Bose gas with 
attractive boundary conditions and the inhomogeneous condensate localization. 

Finally we repeat that we analyzed only the problem taking the thermodynamic limits 
in the sense of homothetycally increasing cubes, with the consequence that the condensate 
is situated anisotropically in the direction of the corners of these cubes and is "localized 
at infinity". To prove this we tune the homothety point position at the logarithmic (in 
the units of the cube size) distance from the cube boundary. If instead one looks for the 
limit of spherical containers, this anisotropy in the positioning of the condensate should 
disappear. Does one expect spontaneous spherical symmetry breaking of the equilibrium 
states in this case? 
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